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Abstract. The lower bound on the period of uniform rotation of neutron stars 
(NS) with causal equation of state (EOS), P m i n , is shown to be determined, with 
a high precision, by the upper bound on the surface rcdshift of static neutron star 
models, z max . We obtain z max = 0.851, which turns out to be about 5% lower than 
upper bound obtained by Lindblom (1984). Using then the empirical formula for the 
maximum rotation frequency of uniform rotation of neutron star models with realistic 
causal EOS, derived by Lasota et al. (1996), we establish an approximate relation 
between z max and P m j n , in which P m i, n is proportional to the maximum mass of static 
NS, M™ ax . In this way we reproduce, within 2%, the exact formula of Koranda et 
al. (1997) based on extensive exact numerical calculations. Replacing A/™ 3 " by the 
maximum measured NS mass, 1.442 Mq, we get a lower bound P m i n ~ 0.288 ms. 



1. Introduction 

The lower limit on the period of a uniformly rotating neutron star depends on the - 
largely unknown - equation of state (EOS) of dense matter above the nuclear density. 
Therefore, an uncertainty in the high density EOS implies a large uncertainty in the 
minimum period of uniform rotation, P m ; n (see, e.g. [|], Q, [§, §)■ Hence, it is, 
therefore, of interest to find a lower limit on P m i n , that is independent of the EOS. 
This limit results from the condition of causality, combined with the requirement that 
EOS yields neutron stars with masses compatible with observed ones (currently the 
highest accurately measured neutron star mass is M^| x = 1.442 M Q Q). It will be 
hereafter referred to as P^hr 

The first calculation of P^in was done by Glendenning 0, who found the value 
of 0.33 ms. Glendenning [Q, however, used a rather imprecise empirical formula, to 
calculate a lowest P m i n by using the parameters (mass and radius) of the maximum 
mass configurations of a family of non-rotating neutron star models. His result, 
therefore, should be considered only as an estimate of P^in- Recently, Koranda et 
al. J^] extracted the value of P^in fr° m extensive exact calculations of uniformly 
rotating neutron star models. They have shown, that the method of Glendenning M 
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overestimated the value of P^ n by 6%. The result of Koranda et al. || calculations 
can be summarized in a formula 

/Vf max 

= 0.196 -f - ms , (1) 

which combined with measured mass of PSR B1913+16 yields today's lower bound 
for Pj~^ = 0.282 ms. This absolute bound on the minimum period was obtained for 
the "causality limit (CL) EOS" p = (p — po)c 2 , which yields neutron star models of 
the surface density po and is maximally stiff (dp/dp = c 2 ) everywhere within the star; 
it does not depend on the value of po- In the present contribution we show that Eq. 
(1) can be reproduced using an empirical formula for P m i n derived for realistic causal 
EOS, combined with an upper bound on the relativistic surface redshift for static 
neutron stars with causal EOS. 



2. Relation between x s and P m ; n 

Numerical results for maximum frequency of stable uniform rotation of neutron stars, 
obtained by Salgado et al. ||, Q for a broad set of realistic, causal EOS of dense 
matter, can be reproduced within better than 5% by an empirical formula Q 

( GM \ 

simple empirical formula : S] maj{ ~ 0.673 ( — 5^ ) s _1 , (2) 

V / 

where M s and R s are, respectively, mass and radius of static configuration with 
maximum allowable mass. Numerical prefactor in Eq.(^|) is independent of the EOS 
from the set considered in ^) , B. 

Defining a relativistic compactness parameter 

_ 2GM B 

related to the surface redshift of static configuration with maximum allowable mass 

by 

Zs = (i-x s y 1/2 -1 , (4) 

we can transform Eq. (^) into an approximate empirical formula, expressing the 
minimum period of rotation of neutron star, P m j n = 27r/O max , in terms of M s and x s , 

simple empirical formula => P min **™**L ms . (5) 

xj M © 

Clearly, for a given maximum mass of static configurations, M s , the lower bound on 
Pmin will be obtained by replacing x s by its upper bound, x S!max . 

The formula for P m ; n can be improved by using more precise empirical formula 
for f2 max . Such a formula was constructed by Lasota et al. fIo|| . As shown in fl0|| , 
numerical results of || , Q for the maximum frequency of uniform stable rotation can 
be reproduced (within better than 2%), for a broad set of realistic causal EOS of dense 
matter, by an improved empirical formula 

<GM ^ 1/2 



improved empirical formula : f2 max C(x s ) f ) s 1 , (6) 
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where C(x s ) is a universal (i.e. independent of the EOS) function of the compactness 
parameter x s , 

C(x s ) = 0.468 + 0.378a; s . (7) 
Combining Eq. (g) and Eq. (|) we get 

A •• , f , p 8-754 x IP' 2 Af s 
improved empirical formula =>• -r m in — 57^ — ms. (8) 

C(x s )xl' 2 M 

Both "empirical formulae" for P m ; n , Eq.(||) and Eq. (|8]) , imply, that at given maximum 
mass of a spherical configuration, the minimum rotation period is obtained for the 
maximum value of x s . At fixed x s , the value of P mm is proportional to M s . Neutron 
stars for which masses have been measured, rotate so slowly that their structure can 
be very well approximated by that of a spherical star. Observations impose thus a 
condition M s > M™f. 



3. Upper bounds on surface redshift 



To minimize P m j n for given M^ x , we have to look for an EOS which yields maximum 
x s at M s = M^ x . It is well known, that if one relaxes the condition of causality, 
the absolute upper bound on x s for stable neutron star models is reached for an 
incompressible fluid (i.e., p = const.) EOS; the value of x s is then independent of M s 
and equal 8/9; it corresponds to an upper bound on z s , equal 2 (see, e.g., ]ll|). 

It is therefore rather natural to expect that in order to maximize z s under the 
condition of causality, one has to maximize the sound velocity throughout the star. 
Together with condition of density continuity in the stellar interior this points out at 
the CL EOS, p = (p — po)c 2 , as to that which yields "maximally compact neutron 
stars"; introducing density discontinuities does not increase the value of z s , see |l2| . 
Note, that the value of z s for CL EOS does not depend on po (and therefore is M s - 
independent). It represents an absolute upper bound on z s for causal EOS, z maK . 
Our numerical calculation gives z s (CL EOS) = z max = 0.8509. This corresponds to 
an absolute upper bound on the compactness parameter of neutron star models with 
causal EOS, x max = 0.7081. 

Let us consider the effect of the presence of a crust (more generally, of an envelope 
of normal neutron star matter). For a given EOS of the normal envelope, the relevant 
(small) parameter is the ratio pb / PbC 2 , where ph and p\, are, respectively, pressure and 
mass density at the bottom of the crust |l3| . The case of pb = corresponds to stellar 
models with no normal crust. Numerical calculations show, that adding a crust onto 
a CL EOS core implies an increase of R s , which is linear in Pb/PbC 2 ; for a solid crust 
we have typically Pb/pbC 2 ~ 10~ 2 . The change (increase) in M s is negligibly small; it 
turns out to be quadratic in ph/pbC 2 . This implies, that the decrease of x Sjmax , and of 
the maximum surface redshift z s ,max, due to the presence of a crust, is proportional to 
Pb/pbC 2 . This is consistent with Table 1 of Lindblom [[13"| . However, the extrapolation 
of his results to pb = yields z max = 0.891, which is nearly 5% higher than our value 
of ;z m ax ! This might reflect a lack of precision of the variational method used by 
Lindblom |TJ|, which led to an overestimate of the value of z ma x- It should be stressed 
that while a precise determination of M max = M s for static neutron star models is 
rather easy, determination of the precise value of the radius of the maximum mass 
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configuration, i? s , (with the same relative precision as M s ) and consequently of the 
value of x s (with, say, four significant digits), is much more difficult and requires a 
rather high precision of numerical integration of the TOV equations. 

4. Lower bound on P m i n 

We restrict ourselves to the case of the absolute lower bound on P m i n , obtained for 
neutron star models with no crust. Inserting the value of £ s ,max into Eq. (|^) we get 

simple empirical formula ==> P^± n ~ 0.22 — ms . (9) 

M 

The current lower bound on P, resulting from the above equation, is 0.32 ms, which 
should be compared with the lower bound of 0.282 ms, obtained in extensive exact 
numerical calculations by Koranda et al. ||. 

Using the improved empirical formula, Eq. (g), we get 

M max 

improved empirical formula P^t n — 0.1997 — ms ^ ^0) 

M Q 

which implies a current lower bound of 0.288 ms, only 2% higher than the exact result 
of Koranda et al. [||. 

The formula (FLO) deserves an additional comment. In numerical calculations, of a 
family of stable uniformly rotating stellar models, for a given EOS of dense matter, 
one has to distinguish between the rotating configuration of maximum mass, which 
corresponds to the rotation frequency r2^/ max (EOS), and the maximally rotating one, 
which rotates at f2 max (EOS) [gj, Notice, that determination of a maximum 

mass rotating configuration (and therefore of ^A/ max ) is a much simpler task than the 
calculation of exact value of f2 max , which is time consuming and very demanding as far 
as the precision of numerical calculations is concerned. Usually, both configurations 
are very close to each other, and il max is typically only 1-2% higher than fiM max ; such 
a small difference is within the typical precision of the empirical formulae for f2 max . 
Actually, the formula for C(x s ), Eq. (||), was fitted to the values of f^/ max (EOS) 
calculated in ||, §]. Therefore, Eq.([lO|) should in principle be used to evaluate the 
causal lower bound to -P m in,M max ; it actually reproduces, within 0.2%, the exact 
formula for this quantity, obtained by Koranda et al. || [see their Eq. (8)]. 

General experience shows that - in contrast to interpolation - extrapolation is 
a risky procedure. The fact that in our case extrapolation of an empirical formula 
yields - within 2% - the value of P m ; n of Koranda et al. || (and reproduces their 
value of -Pmin.Afmax), proves the usefulness of compact "empirical expressions" which 
might summarize, in a quantitative way, a relevant content of extensive numerical 
calculations of uniformly rotating neutron star models. 
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